Phonon Softening in Superconducting Diamond 
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We observe strong softening of optical phonon modes in superconducting (T c = 4.2 K) boron- 
doped diamond near the Brillouin zone center using inelastic x-ray scattering from a CVD-grown 
highly oriented sample. The magnitude of the softening, and our observation that it becomes 
stronger approaching zone center, supports theoretical models suggesting a phonon-mediated pairing 
mechanism via coupling of optical phonon modes to Fermi surfaces around the zone center. The 
electron-phonon coupling parameter is determined as approximately A = 0.33. 



Pure diamond is a wide bandgap insulator. When it 
is boron-doped beyond the metal-to-insulator transition 
at riB ~ 3 • 10 20 cm -3 , it shows superconductivity with 
transition temperatures T c that are remarkably high for 
an impurity doped material with low carrier density. Su- 
perconductivity was reported both in samples grown by 
high-pressure-high-temperature synthesis with T c = 4 K 
at a B-doping level n# ~ 5 • 10 21 cm -3 [1 and in sam- 
ples grown by /i-wave assisted chemical vapor deposition 
(CVD) [2j [3]. T c depends on the doping level tlb and 
increases up to T° n w 11 K (onset) at the highest doping 
reported so far [4]. 

Conventional theory of phonon-mediated superconduc- 
tivity suggests that an ultrahard material like diamond 
can have a high T c if the high frequency modes have 
a strong enough electron-phonon coupling (EPC). The 
bond-stretching optical modes of diamond around c^o = 
160 meV naturally couple to the covalent bonding states 
of the diamond valence band. Various ab-initio mod- 
els have predicted the EPC coupling constant A. In the 
typical doping range between n = 2.5 — 3 at%, values be- 
tween between A = 0.27 (virtual crystal approximation 
VCA [5]) and A = 0.39 (supercell [6]) are calculated and 
lead to predicted T c = 0.2 — 4.4 K that correspond well to 
those observed in real samples. The nature of the modes 
that couple to the free charge carriers is, however, diffi- 
cult to assess by theory and few experimental data are 
available. VCA naturally focuses on the optical modes 
of the diamond lattice that couple due to their bond- 
stretching nature [5j [7] . The A calculated in a supercell 
also includes localized B-modes, where the coupling is en- 



hanced due to the enhanced valence band wave functions 
at the dopant sites [8]. 

We have studied the lattice dynamics in a sample of 
highly B-doped superconducting diamond. For the opti- 
cal phonon branches, we observe a softening that allows 
the determination of the momentum-dependent EPC pa- 
rameter X(q). Similar to inelastic neutron scattering, the 
technique of inelastic x-ray scattering (IXS) probes the 
momentum-dependent frequency of the phonon excita- 
tions. A strong softening of the optical phonon branches 
is the fingerprint of a strong EPC. Moreover, the ob- 
servation of the momentum-dependence of the softening 
provides information on the phase-space of the coupling. 
This complements early indications of phonon soften- 
ing from Raman scattering experiments that probe only 
zero momentum [9]. In the case of high B-doping, Ra- 
man scattering shows a strong signal from both lattice as 
well as electronic excitations that could not be disentan- 
gled [10] and thus no quantitative data on the EPC have 
been gained from the Raman spectra. A recent study of 
far-infrared and THz-radiation absorption spectra shows 
one peak of EPC in the region of the optical phonons 
as well as contributions in a region of lower frequency 
modes [TT] . 

The electronic structure in the metallic regime is de- 
scribed either by a shift of the Fermi level into the di- 
amond valence band (degenerate metal) or by a purely 
B-related impurity band or a mixture of both. Experi- 
mentally, a Fermi level in the diamond valence band was 
observed by angle-resolved photoemission spectroscopy 
(ARPES) [12 , but evidence for a persistence of an impu- 
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FIG. 1: (a) Schematic view of the ~ 100/xm thick B-doped 
sample B. (b) Temperature-dependent resistance curve. 
Above T c the resistance is almost flat with a room tempera- 
ture resistivity of ~15 Hem. (c) X-ray rocking curves of the 
(004) reflection of sample B (o) and undoped sample N (a). 



rity band close to the top of the valence band maximum 
was reported by optical spectroscopy [13]. The former 
scenario would lead to a well-defined "spheroid" Fermi 
surface (FS) around the T-point. This FS has signifi- 
cant anisotropy due to the covalent bonding states and 
consists of three concentric sheets. The latter scenario 
should not lead to a well-defined FS if the dopant atoms 
occupy random positions in the lattice. Calculations us- 
ing the coherent phase approximation suggest the forma- 
tion of a spheroid FS around T but a smearing of the 
Fermi wave vector due to impurity scattering [13] . 

In the presence of a spheroid FS a strong phonon soft- 
ening is expected at small momentum g, where energy- 
conserving phonon scattering is possible within the FS. 
The softening is gradually reduced to a constant value at 
increased phonon momentum around the Brillouin zone 
boundary. This gentle kind of a Kohn anomaly is rem- 
iniscent of the EPC in MgB2, where the Cooper pairs 
form due to a coupling of the electronic cr-band to the 
E2 g optical phonons [5] [15j [16] . 

The disorder in B-doped diamond has only a little di- 
rect effect on the lattice dynamics, because the mass dif- 
ference of B- and C-atoms is small. Like the natural 
isotope mixture of C-atoms in diamond it would have a 
significant effect only at doping levels much higher than 
a few atomic percent [17] . In our data we found evidence 
for a slight disorder-induced change of the Eigenvectors 
of acoustic modes [18] . but no evidence for an EPC of 



acoustic or localized modes was found. In this letter we 
discuss the coupling of the optical phonons, where a di- 
rect comparison between pure and doped diamond leads 
to a quantitative determination of the mode-dependent 
coupling parameter A. 

The samples were grown by heteroepitaxial CVD depo- 
sition on a Si(001) substrate pre-treated by SiC to yield 
a highly ordered (OOl)-oriented ensemble of coalesced di- 
amond grains [19] . A layer of 100 /im B-doped diamond 
at a doping level of ub — 3.8 • 10 21 cm -3 w 2 at% was 
deposited in ~100 hours (sample B). The grains are ori- 
ented within 1° as seen in the x-ray rocking curve in 
Fig. [1J. The superconducting transition is rather broad 
around T c = 4.2 K (see Fig. [TJd). The silicon sub- 
strate was removed by chemical etching. A second sam- 
ple of 100 /im of more lightly nitrogen-doped diamond 
(ti 7v = 4 • 10 18 cm -3 ) was grown by the same technique 
(sample N). Nitrogen is introduced to enhance the CVD 
growth. It does not lead to an electronic doping due to 
its deep impurity level. Our experiment confirmed the 
phonon dispersion of sample N to be identical to pure 
diamond [20] and we will refer to this reference sample 
as undoped. The doping levels were determined by sec- 
ondary ion mass spectroscopy (SIMS). 

The experiment was performed at the IXS station of 
BL35XU [21] at the SPring-8 synchrotron light source. 
X-rays of hv = 15.82 keV were selected with energy res- 
olution of 6.4 meV by use of the Si (888) reflection in the 
backscattering monochromator and in the diced spheri- 
cal analyzers. For each of the two samples 8 positions 
of the diffractometer arm were measured, which results 
in 8 spectra precisely on the high symmetry lines T-L 
and T-X. Additional 21 spectra very close to the high 
symmetry line at identical positions for undoped and B- 
doped diamond were also analyzed. These spectra show 
only longitudinal phonons (LO) due to the measurement 
geometry. In addition, a spectrum was acquired at the 
L-point for the transverse phonons (TA and TO). The 
energy scale was calibrated for each analyzer using the 
elastic line position and the optical phonon peak from a 
pure single crystal diamond sample close to the T-point 
at Q = (1.1 1.1 1.1) [2D] . All data were acquired at room 
temperature because the EPC as such has very weak tem- 
perature dependence and can be studied much more con- 
veniently at room temperature than at lower tempera- 
ture. We remark at this point that even at temperatures 
below T c no change of the EPC is expected for the opti- 
cal modes that have frequencies c^o much larger than the 
superconducting gap c^o ^> 2 A w 1 meV [TT| 122]. 

A selection of measured spectra in the region of the 
optical phonons is shown in Fig. [2] For small momenta 
close to the T-point, the peaks are strongly shifted with 
respect to undoped diamond (extrapolated —8 meV). At 
the Brillouin zone boundary the shift is smaller (—2 meV) 
both for the X-point and the L-point and for LO as well as 
for TO phonons. At the zone boundary the peak shape is 
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FIG. 2: IXS spectra from undoped sample N (v) and sample B (o) for momentum transfer along the [001]-direction T-X (a and 
b) and the [111] direction T-L (c and d) in a geometry where only longitudinal phonons contribute, as well as at Q = (0.5 1.5 0.5) 
(e) where mostly transverse phonons contribute. The total momentum transfer Q is indicated for each pair of spectra in relative 
lattice units. The peak position and width at half height is indicated. 



given by the resolution function of the IXS spectrometer. 
Near the zone center the peaks from sample N are still 
resolution-limited, while the peaks from sample B are 
very broad and asymmetric. This asymmetric peak shape 
is incompatible with a life-time reduction of the phonons 
in a damped harmonic oscillator theory. More probably 
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FIG. 3: (a) Dispersion of the optical phonons in sample N (v) 
and B (o) as determined by the peak-positions as shown in 
Fig. [2] In the left-hand side the T-L and in the right-hand 
side the T-X direction is shown. Most data are from the LO 
phonons, but in the left-hand side at L the TO phonons are 
also indicated, (b) Softening curve: peak shift for sample 
B with respect to sample N as a function of momentum, (c) 
Momentum-dependent EPC parameter X(q) as determined by 
Eq. (§. 



it is given by a certain degree of in-homogeneity of the B- 
content across the sample that is also evident in the broad 
superconducting transition shown in Fig. [TJd. The same 
peak positions and asymmetric shapes were recently also 
observed in a sample grown by homo-epitaxial C VD [23J . 

The dispersion shown in Fig.[3^ is determined from the 
peaks as the arithmetic middle between the two flanks 
at half height as indicated in Fig. [2] The difference of 
the dispersions for samples N and B is shown in Fig ^jp. 
Towards T, the peak shift approaches an extrapolated 
value of 8 meV. This value is larger than the one naively 
derived from the shift in Raman spectra (< 4 meV), 
where the one-phonon line remains barely visible above 
the background of the lower-energy Fano- feature. The 
momentum-dependence of the softening is compatible 
with a gentle Kohn Anomaly in the spheroid FS of a 
degenerate metal, as the two directions and also longi- 
tudinal and transverse phonons soften in the same way. 
Thus, we find that the phonon softening is purely depen- 
dent on the absolute value of the momentum q for the 
two high symmetry direction we investigated. 

To determine the ^-dependent EPC parameter A(g), 
the first approach would be to consider the life-time 
broadening of the phonons and use Allen's formula that 
connects X(q) to the peak width 7(g). In diamond, this 
approach fails because of the in-homogeneity induced 
broadening of the LO peaks near T and their asymmet- 
ric line-shape. Instead, we can determine X(q) from the 
softening from Allen's analysis [24] as 

where E c is a cut-off energy identified with the Fermi en- 
ergy E c = Ep, N(0) is the density of states at Ep (in 
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states/eV/cell) and uo(q) is the measured dispersion for 
samples N and B. We take N(Q) and Ep from VCA calcu- 
lations for our doping level (N(0) = 0.12 states/eV/cell 
and E F = 0.62 eV) [7 . The product N(0)E F is pro- 
portional to the number of holes per unit cell (in a 
parabolic band in diamond N(0)E F = 3n), so that dif- 
ferent models of the band dispersion still give only a 
small error in the determination of A (q). The relevant 
EPC parameter is the average over the Brillouin zone 
A = 1/Vbz J bz X(q)d 3 q. Assuming isotropic softening 
we can extrapolate X(q) from either of the two measured 
momentum directions independently and determine the 
Brillouin zone average as ArL = 0.39 and Arx = 0.27, re- 
spectively. The difference gives an approximate measure 
of the experimental uncertainty of this analysis. The final 
value is the average between the two: A = 0.33 ± 0.06. 
The systematic uncertainty due to the theoretical ap- 
proximations is probably of the same magnitude (±20%). 

The experimental EPC parmeter lies in the range of 
theoretical predictions. When inserted into the McMillan 
relation T c = (u /1.2) exp [-1/ (A/(l + A) - /i*)] using 
the optical phonon frequency uoq = 160 meV and a con- 
ventional /i* = 0.12, a T c = 0.63 K is derived (T c = 3 K at 
the higher error limit of A), which is significantly smaller 
than the observed T c = 4.2 K. In a direct comparison the 
measured softening is less strong than the predictions by 
VCA calculations (~30 meV at T and ~6 meV at the 
L-point [25]). Supercell calculations that provide a more 
realistic incorporation of the B-atoms into the diamond 
lattice suffer from ordered-boron artifacts, but the soften- 
ing reported a a higher doping level of 6.25 at% (12.5 meV 
at T [6 ) is in better agreement with our observations. In 
addition to the optical phonons, these calculations also 
include coupling to localized B-modes. Thus, although 
the observed coupling of the optical phonons is strong, 
there is room for other modes to provide additional EPC. 

From our data the determination of A is only possible 
for those modes that have a direct relation between the 
doped and the undoped case, i.e. the six phonon branches 
of the diamond lattice. No coupling was observed for 
the acoustic modes [18]. Also, we have not observed a 
clear signature from B-derived, localized modes. In the 
IXS spectra we have, however, observed finite intensity 
almost throughout the spectral range that is probably 
related to impurity modes. Our data are thus consistent 
with a recent optical spectroscopy experiment that found 
significant contributions to the EPC both in the high- 
frequency spectral range above uo = 140 meV as well as 
in the lower range around uo = 70 meV [11 . The former 
region corresponds well to the optical phonon coupling 
reported in this letter. The latter region is also the range 
of B-B pair vibrations that lead to a very large Raman 
scattering signal but do not contribute to the EPC [26] . 
Contributions from single boron localized modes are also 
predicted around 70 meV [6 , 26 , matching the analysis of 
Ref. [11 , but so far unaccounted by direct spectroscopy. 



In conclusion, we have measured the softening of the 
optical phonons in diamond. We found a strong soft- 
ening with a characteristic momentum dependence both 
along T-L and T-X. This momentum dependence is com- 
patible with coupling to the holes in a spheroid FS as 
predicted by theory for a Fermi level in the valence band 
of diamond. We have determined the momentum depen- 
dent EPC parameter X(q) for the optical phonons and 
its average over the Brillouin zone is the mode-specific 
coupling parameter A = 0.33 ± 0.06. Our data thus 
strongly support an Eliashberg model of the supercon- 
ductivity in B-doped diamond with coupling through the 
bond-stretching optical modes that are particularly high 
in frequency for the ultrahard diamond material. 
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